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I ' An alternative approach to symmetry restoration within Energy Density Functional, the 

^--s I Symmetry-Conserving EDF is discussed. In this approach, the energy is directly written 

^sj . in terms of the degrees of freedom encoded in the one-, two-... body density matrices of 

the state having good quantum numbers. The SC-EDF framework is illustrated within 

i2 ' Projection After and Before variation applied to particle number restoration. 
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§1. Introduction 

The most popular method to describe pairing correlations in many-body systems 
is based on the Bardeen-Cooper-SchrieflFer (BCS) or Hartree-Fock-Bogoliubov (HFB) 
approaches where one uses quasi-particle states that explicitly break the U(l) sym- 
metry associated with particle number conservation. This approach is extensively 
y ■ used in Energy Density Functional (EDF) theory. It is, however, known that the 

rH . BCS or HFB description of pairing can be greatly improved, especially in the regime 

of weak coupling and/or small particle number, if the particle number is properly 
restored.!^ Such a restoration is usually made using configuration mixing technique 
^ . corresponding to a so-called Multi- Reference EDF (MR-EDF). 

Jl^ I BCS, HFB and symmetry restoration are well-defined techniques when using a 

Qs ■ Hamiltonian operator. Additional care has to be taken when using a more general 

rn ■ EDF eventually based on effective interactions. The validity and range of applica- 

(^ \ bility of configuration mixing within EDF methods has been debated recently.'^ ™ 

^^ ■ In particular, (i) the necessity to regularize the MR-EDF,'^ ''^'"^'' (ii) the density to 

be used in density dependent interactionjl^ and (iii) the need to clarify the notion 
of breaking/restoration of symmetries in density functional theory^ have been dis- 
cussed. 

p\ . Here, the restoration of particle number within Energy Density Functional the- 

^ I ory is re-analyzed. We show that the MR-EDF can, under certain conditions, be 

interpreted as a functional theory of the one-body and two-body density matrices of 
the projected product state. Starting from this observation, we propose a new density 
functional approach, which we call Symmetry-Conserving EDF (SC-EDF) where the 
breaking and restoration of symmetries are accounted for in a single step-l^ This ap- 
proach has the advantage to provide a suitable framework that solves the difficulties 
(i-iii) and to be directly applicable to nuclei using most popular effective interaction. 
Applications to light and medium mass nuclei of Projection After Variation (PAV) 
and Variation After Projection (VAP) will be presented. 
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§2. The nuclear Energy Density Functional 

In Density Functional Tlieory (DFT) or Energy Density Functional (EDF) ap- 
plied to nuclei, one generally starts from a set of trial quasi-particle states, and writes 
the energy in terms of underlying densities. For instance, for the effective Skyrme 
or Gogny interactions, the EDF can be written as: 

£ [^o] =£[p,K,K*] = Y^ tup,, + - ^ v^i-ijPiiPjj + lYl "^m^^Sj ' (2-1) 

where p and k are the normal and anomalous density of the state ^o respectively, 
whereas v^^ and v'^'^ are the effective, eventually density dependent, kernels in the 
mean-field and pairing channel, respectively. 

At present, EDF methods are the only theory that provides a unified framework 
to describe nuclear structure, reaction and thermodynamics. One of the keys for this 
success is the use of trial states that do explicitly break some of the symmetries of the 
underlying many body Hamiltonian: translational invariance, parity, U(l) symmetry 
associated to particle number, ... Thanks to symmetry breaking, correlations like 
pairing can be included with a rather simple functional. 

The EDF based on a single trial state (Single-Reference [SR] approach) can 
provide a suitable approach for many aspects of nuclear physics. In applications to 
nuclear spectroscopy however, it is often necessary to have states with good quantum 
numbers. 

The natural way to tackle this problem is to restore symmetries that have been 
broken at the SR level. Restoration of broken symmetries within nuclear EDF is 
strongly guided from the configuration mixing (Generator Coordinate Method) ap- 
proach generally used in the Hamiltonian case (see discussion in'l''). This leads to 
the so-called Multi-Reference concept (MR-EDF). Unless a Hamiltonian is used as 
a starting point to design the EDF and no approximation is made on the exchange 
term, the justification of GCM like approaches is unclear. Recently, many practi- 
cal and formal difficulties have pointed out some weakness of MR-EDF as they are 
currently formulated.'^ ™'Bli -13. Eli in particular, after restoration, the new functional 
may contains spurious components that do not properly behave under the symme- 
try group transformations. This has been clearly illustrated for the case of Particle 
Number Restoration (PNR) case^'l^ and further discussed in ref.® for the case of 
rotational invariance. 

§3. The Symmetry-Conserving EDF concept 

We recently have shown that restoration of broken symmetries can eventually be 
made using a slightly different strategy and introducing the concept of Symmetry- 
Conserving EDF. This strategy is presented here. Let use consider an EDF given by 
Eq. (|2-ip and assume that a given symmetry <S (spherical symmetry, particle num- 
ber, ...) is broken. Introducing the transformation associated with this symmetry, 
denoted by R{Q). This operator can for instance represent a rotation in 3D space 
for the case of broken spherical symmetry or rotation in gauge space for the case of 
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the U{1) symmetry associated with particle number. Starting from (pQ, the energy 
of any rotated state, denoted by \'P{f2)) = i?(J7)|^o) is generally degenerated with 
the original energy £[<l>{f])] = £[(Pq]. To restore broken symmetries, a new trial state 
is generally introduced that sums up the different orientation, i.e. 

\^s) = Ps\.Pq) = f dnG{n)R{n)\<Po) = [ dnG{n)\'P{n)) (s-i) 

Jn J Q 

where G{Q) are weight factors that depend on the symmetry to be restored, and Ps 
is the projector on good quantum numbers associated with the symmetry. 

The standard approach, guided by the GCM technique, to restore symmetries 
within EDF is to directly write the energy as 

£s[M = I dQ£ rpO^,KO^,K^o*l A/'5(0,^) . (3-2) 

where 7V5(0, i?) is a factors depending on G{f2), while p^ , k^ , ... are the transition 
densities associated with the couple of states {^0;^(^)}- Defining the energy as in 
(|3-2|) for an EDF with different vertices v^^ and v'^'^, that in addition might not be 
completely antisymmetric and contain elaborate density dependencies, will lead to 
spurious terms that might eventually spoil MR-EDF applications. 

Recently, we proposed to write the energy directly in terms of the degrees of 
freedom associated with the projected state,'^ i.e. 

£sWs]=£s[p^,R^,...]. (3-3) 

where p^ , R'^ , ... denote, respectively, the one-, two-, ... body densities of a sym- 
metry restored state. Such an approach called hereafter symmetry conserving is 
motivated by the following reasons: 

• When a Hamiltonian is used equation (|3-3|) and (|3-2|) are strictly equivalent.'^ 

• For functionals that can be regularized, it has been shown recently that, in the 
case of PNR, the regularized MR-EDF takes a form very close to (j3-3p ^ 

• Since the building block of the functional theory are observables of the projected 
state, they automatically include the proper properties under the transforma- 
tion of the restored symmetry and so do the energy after projection. 

• The SC-EDF can be regarded as an extension of the SR-EDF onto a different, 
eventually larger Hilbert space of trial states. 

Illustrations of the SC-EDF application are given below for the PNR case. 

§4. Applications 

Starting from a general discussion on MR-EDF applied to particle number 
restoration, it has been shown that the SR-EDF can be conveniently extended to 
include the conservation of particle number. In that case, the energy reads 



U^N] = Y.iunf^ + l E ^.^i^. + i E ^m^. 
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Fig. 1. (Color online) Left: Particle-number restored deformation energy curve of ^*0 as a function 
of quadrupole deformation Pi and calculated within standard MR-EDF technique using Pro- 
jection After Variation (PAV) with the parameterization SIII of the Skyrme EDF and a delta 
pairing interaction without (blue dashed curve) and with regularization (black solid curve) . The 
red filled circles correspond to the result obtained using directly Eq. (|4-l|l (see text). To compare 
with previous work,^' the Coulomb exchange contribution has been subtracted from the energy. 
Right: Same for the Sly4 parameterization in the mean-field channel. In that case, only two 
curves can be shown since the MR-EDF cannot be regularized. To emphasize the steps and 
divergences, the MR-EDF calculations were performed with many more discretization points of 
the gauge-space integrals than usually done, cf ref.'^ for further details. 
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where n^ denotes the occupation numbers and R^ the two-body density matrix of 
the state with particle number N , denoted by \^n)- 

This SC-EDF has been used to perform Projection After Variation (PAVJ^J 
similarly to what is currently done in MR-EDF (see figure [T]). When regularization 
of MR-EDF is possible, the energy obtained using Eq. ()4-ip is strictly equivalent 
to the result of the regularized EDF (left of figure [1]). Contrary to MR-EDF, the 
SC-EDF can also be applied to functional containing a non-integer powers of the 
density (right side of [1]) leading to a smooth behavior of the potential energy curves. 

More recently, the SC-EDF has been used to perform Variation After Projection 
(VAP). In that case, the energy is minimized directly by making the variation of the 
projected state degrees of freedom. In figure[2]an illustration of the energy gain in the 
pairing channel is shown for the krypton isotopic chain. One of the main advantages 
of VAP is to give also non-zero pairing energy for closed-shell nuclei. For instance, 
the ^^Kr is at a neutron sub-shell closure N = 40. For the pairing functional used 
here, HFB or BCS approaches lead to zero pairing energy in the neutron channel, 
whereas the VAP gives a pairing energy of 4.1 MeV. Note that this energy gain is 
partially compensated by the rearrangement of the mean-field, leading to a net gain 
of energy of 1 MeV. 



Application of symmetry breaking and its restoration in finite nuclei 




35 40 45 

Neutron number 



50 



Fig. 2. (Color online) Illustration of the (proton+neutron) pairing energy obtained after PNR-VAP 
(solid curve). The dashed curve correspond to the original SR-EDF result. 

§5. Summary 

Motivated by recent difficulties observed in MR-EDF, we have recently intro- 
duced a new symmetry-conserved EDF framework. Application to particle number 
restoration are very encouraging. In particular, this approach can be applied to func- 
tionals with rather general dependence on the density, unlike standard regularized 
MR-EDF. In addition, its application to VAP, illustrates its feasibility. The main 
difficulty is now to provide functional forms that can be used for other symmetry 
restoration like deformation. 
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